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摘 要 : 为 揭示 银川 平原 降水 稳定 同位 素 的 时 间 尺 度 效应 及 水 汽 来 源 ,采用 后 向 轨迹 模 型 聚 类 分 析 
法 潜在 源 贡 献 分 析 法 (PSCF) 和 浓度 权重 轨迹 分 析 法 (CWT) 解 析 水 汽 来 源 及 洪 在 蒸发 源 区 。 结 果 
表明 :(1) 银川 平原 降水 毛 氧 稳定 同位 素 组 成 存在 明显 的 季节 性 变化 , 冬 半 年 降水 氨 氧 稳定 同位 素 
组 成 (-38.6%o+51.6%o 和 -4.5%o+5.2%o) 明 显 偏 正 于 夏 半年 (-40.9%o+17.7%o 和 -5.7%o+3.0%o) ; 不 同 
时 段 下 ,局 地 大 气 降水 线 (LMWL) 具 有 显著 差异 , 夏 半年 LMWL 的 斜率、 截 距 (5.43、-9.71) 低 于 冬 半 
年 的 斜率 、 截 距 (9.10、5.08) 和 银川 平原 全 年 LMWL 的 斜率 、 截 距 (6.79、-2.79)。(2) 在 全 年 、 夏 半年 
和 冬 半 年 ,降水 中 8*0 的 气温 效应 显著 ,气温 效应 分 别 为 (0.473+0.210)%o*%C ', (0.258 + 
0.037 )%o + °C 和 (0.211+0.031)%o.%C ,降水 量 效应 不 显著 。(3) 在 夏 半年 和 冬 半 年 ,研究 区 降水 的 
水 汽 来 源 主要 为 西风 水 汽 , 同 时 受 局 地 蒸发 水 汽 的 影响 。 夏 半年 ,水 汽 蒸发 源 区 主要 分 布 在 研究 


区 周边 地 区 及 西部 、 


南部 和 南部 ; 冬 半年 主要 分 布 在 研究 区 周边 地 区 和 和 西部。 研究 结果 可 为 干 


旱地 区 利用 稳定 同位 素 技 术 鉴 别 水 汽 来 源 方法 的 选择 提供 理论 依据 。 
X 键 词 : 大 气 降水 ; 稳定 同位 素 效应 ; 水 汽 来 源 ; 潜在 蒸发 源 区 ; 不 同时 间 尺 度 ; 银川 平原 


文章 编号 : 


银川 平原 地 处 我 国 西北 干旱 区 ,水 资源 是 影 
响 和 制约 该 地 区 高 质量 发 展 和 生态 建设 的 关键 因 
素 沁 。 降 水 是 水 循环 过 程 中 的 一 个 重要 环节 ,研究 
降水 稳定 同位 素 的 时 间 尺 度 效应 ,可 有 效 示 踊 区 域 
水 汽 来 源 ,定量 阐明 降水 在 水 循环 过 程 中 的 分 配 和 
转化 2 ,对 进一步 了 解 该 区 域 水 循环 过 程 具有 重要 

降水 氨 氧 稳定 同位 素 组 成 受气 温 .降水 量 和 水 
汽 来 源 等 因素 的 影响 ,随时 间 、 地 区 和 大 气 环流 条 
件 而 发 生变 化 二 。Gui 等 了 对 青藏 高 原 与 干旱 区 过 
渡 区 降水 稳定 同位 素 的 环境 效应 及 时 空格 局 进行 
了 研究 ,发 现 降水 中 8*0 的 气温 效应 为 0.47%o ,在 夏 
季 降 水 量 的 影响 明显 ;Salamalikis 等 利用 浓度 权重 
轨迹 分 析 法 (CWT) 结 合 气 象 及 水 汽 氛 型 余 (d-ex- 
cess) 数 据 , 分 析 了 希腊 帕 特 雷 地 区 的 水 汽 源 区 , 进 


收 稿 日 期 : 2021-03-11; 修订 日 期 : 2021-08-24 
基金 项 目 : 国家 自然 科学 基金 (31760236 ,31460220 ) 资 助 


作者 简介 : 高 阳 (1994- ) , 男 ,硕士 研究 生 , 主要 从 事 同 位 素 水 文 等 方面 


一 步 得 出 水 汽 的 潜在 蒸发 源 区 主要 位 于 地 中 海盆 
地 中 部 地 区 ; 孟 鸿 飞 等 中 利用 TrajStat 软件 中 的 后 
向 轨迹 模型 (HYSPLIT) 和 潜在 源 贡 献 分 析 法 
(PSCF) ,对 黑河 上 游 夏季 降水 的 潜在 蒸发 源 区 进 
行 了 研究 。 现 有 关于 银川 平原 降水 中 稳定 同位 素 
的 研究 30 , 主要 基于 IAEA/WMO 监测 网 中 银川 站 
点 的 月 平均 降水 稳定 同位 素数 据 , 但 受 限 于 监测 站 
点 数量 和 数据 时 效 性 ,还 不 能 完全 了 解 降 水 中 稳定 
同位 素 的 环境 效应 和 水 汽 来 源 。 同 时 ,由 于 研究 区 
地 处 我 国 西北 季风 边缘 区 ,降水 氨 氧 稳定 同位 素 组 
成 特征 不 完全 同 于 西北 其 他 干旱 地 区 ,并 且 水 汽 的 
蒸发 源 区 也 还 不 清楚 。 

本 文 基于 在 银川 平原 银川 和 灵 武 2 个 站 点 实际 
监测 的 日 降水 氢 氧 同位 素 和 相应 气象 数据 ,对 全 
年 、 夏 半年 和 冬 半 年 不 同时 段 下 降水 中 稳定 同位 素 
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的 变化 特征 、 局 地 大 气 降水 线 、 气 温和 降水 量 效应 
进行 了 分 析 和 比较 ,并 探讨 了 该 地 区 降水 的 水 汽 来 
源 及 潜在 蒸发 源 区 ,以 期 为 今后 研究 西北 干旱 区 局 
地 水 循环 过 程 提供 参考 ,同时 提供 可 行 的 水 汽 来 源 
鉴别 方法 。 


1 研究 区 概况 


银川 平原 位 于 黄土 高 原 中 部 西 侧 ( 图 1), 属 于 
温带 大 陆 性 半 干 旱 气 候 , 多 年 平均 气温 为 8.8 C, 
利用 E601 型 蒸 渗 仪 测 得 多 年 平均 蒸发 量 为 1800 
mm, 多 年 平均 降水 量 为 209.7 mm。 降 水 年 内 分 配 
不 均 ,主要 集中 在 6 一 9 月 , 约 占 多 年 平均 降水 量 的 
72.4%( 数 据 来 自 国 家 气象 信息 中 心 - 中 国 气象 数据 
网 提供 的 中 国 地 面 气象 站 逐 小 时 观测 资料 2010 一 
2018 年 数据 )。 


2 数据 与 方法 


2.1 样品 采集 与 测定 

本 文 在 银川 市 宁夏 大 学 贺兰山 校区 (38°30'05”N， 
106°08'34"E ; 764K 1120 m) 和 银川 市 灵 武 市 白 茂 滩 
国家 级 自然 保护 区 内 (38*03'34"N ,106°22'03"E ; Y% 


o 银川 平原 边界 六 
研究 区 
=E 市 行政 边界 


0 
[| 


图 例 aa, A 一 


I 1130 m)2 个 站 点 进行 事件 尺度 的 大 气 降水 样品 
采集 ,降水 样品 明细 见 表 1。 按 照 《 大 气 降水 样品 的 
采集 与 保存 》(CB/T13580.2-1992)5 ,在 降水 日 的 
08:00 Fil 20:00 Xt KE ZK HES 0.1 mm 的 降水 分 别 进 行 1 
次 采样 ,二 者 的 降水 量 之 和 记 为 该 日 的 总 降水 量 。 
在 银川 站 点 ,为 减少 建筑 物 等 其 他 可 能 污染 源 对 降 
水 样品 产生 影响 ,将 2 个 标准 雨量 简 ( 直 径 20 cm 圆 
桶 上 放置 漏斗 ,下 端 与 雨量 筒 相连 ,漏斗 内 放置 乒 
乓 球 防止 蒸发 ) 置 于 楼 项; 在 灵 武 站 点 ,将 标准 雨量 
简 置 于 无 林 空 旷 地 内 。 每 次 降水 结束 后 ,将 收集 到 
的 液态 降水 立即 注入 30 mL 的 聚 乙烯 采样 瓶 中 ,并 
M Parafilm 封口 膜 密封 ;降雪 样品 先 装 入 水 样 袋 中 ， 
密封 袋 口 , 待 其 在 室温 下 完全 融化 后 装 人 30 mL 的 
聚 乙 炳 采样 瓶 中 密封 。 所 有 采集 的 降水 样品 带 回 
实验 室 后 经 0.45 um 的 混合 纤维 滤 膜 过 滤 后 置 于 冰 
箱 中 (4 %C 左 右 ) 冷 藏 保存 待 测 。 

所 有 采集 的 降水 样品 利用 液态 水 同位 素 分 析 
仪 (Los Gatos Research , DLT-100, USA ) 进 行 测定 (5*O0 
FISD 的 测定 精度 分 别 为 +0.1%o 和 +0.4%o)。 测 定 结 
果 的 有 效 位 数 符合 张 自 超 等 ”对 稳定 同位 素 分 析 
结果 的 有 效 位 数 规则 并 以 相对 于 维也纳 标准 平均 
海洋 水 (Vienna standard mean ocean water , V-SMOW ) 


图 1 银川 平原 及 采样 点 分 布 


Fig. 1 Distribution ofYinchuan Plain and samplings points 


表 1 银川 平原 各 站 点 数据 
Tab. 1 Data of the sampling sites in Yinchuan Plain 
采样 点 。 5D 年 均值 /%。 5"0 年 均值 /%。 MAR (d-excess ) 年 均值 /po 平均 气温 /AC ”降水 量 /mm 平均 相对 湿度 /% ”降水 样品 /个 
银川 -29.4 -3.0 -5.7 18.93 101.5 63.85 17 
灵 武 -338 -4.9 5.0 20.84 113.7 70.40 19 
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E HE: 银川 平原 降水 氨 氧 


的 千 分 差 表 示 : 


Ra 
5D( 或 8*0) = sample 


Ry smow 


式 中 :Rw 和 Rvswow 分 别 表示 降水 样品 和 维也纳 标 

准 平均 海洋 水 中 重 轻 同 位 素 丰 度 之 比 (%o)。Dans- 

gaard 7 SC T MÆR (d-excess ) ,可 以 示 踪 降水 的 

水 汽 来 源 , 表 示 为 : 

d- excess = òD - 8 x 8O (2) 

水 汽 中 氢 氧 稳定 同位 素 日 值 由 对 应 的 次 降水 
中 氢 氧 稳定 同位 素 日 值 换算 得 到 ,公式 为 : 

8°O,y =8"0, -10(o 1) (3) 


T Ry -SMOW 


x 1000%0 (1) 


òDpy = 8D, - 10°(a? ,—1) (4) 
式 中 : 8° Opy M SD py 为 水 汽 中 氧 氧 稳定 同位 素 日 降 
水 量 加 权 平 均值 (%o); 8°O, FSD, 为 对 应 的 次 降水 
中 氢 氧 稳定 同位 素 日 降水 量 加权 平 均值 (%o) 5". 
和 a 为 降水 和 水 汽 氧 氧 稳定 同位 素 之 间 的 平衡 分 
馏 系数 ,这 2 个 系数 与 温度 (7,K) 存 在 如 下 关系 “ : 


6 3 
10° hno! = 1.1373) oaiso( 9 ) 2.0667 (5) 


2 


降水 氧 氧 稳定 同位 素 的 加 权 平 均值 计算 公式 
如 下 : 


6 3 
10 Ino? ,= 24.844{ 19 }t6.2as{ 10) +52.612 (6) 


> 5xP 
= Sp 
式 中 :6, 为 降水 量 加 权 平 均值 (%o) ;i 为 降水 事件 ;6; 
为 相对 应 的 稳定 同位 素 值 (%o);P; 为 相对 应 的 降水 
量 (mm)。 
2.2 数据 来 源 及 处 理 
本 文 所 用 气象 数据 为 中 国 气 象 数 据 网 (https:// 
data.cma.cn) 银 川 和 灵 武 站 点 的 逐 小 时 观测 资料 , 数 
据 经 过 质量 控制 ,各 要 素数 据 的 实 有 率 超 过 99.99 ， 
数据 的 正确 率 均 接近 100% ; 气 团 后 向 轨迹 资料 使 
用 美国 国家 环境 预报 中 心 (National Centers for Envi- 
ronmental Prediction ,NCEP) 和 美国 国家 大 气 研 究 中 
心 (National Centers for Atmospheric Research, 
NCAR ) 提 供 的 全 球 再 分 析 资 料及 全 球 资料 同化 系 
统 (GDAS) 气 象 要 素数 据 , 空 间 分 辨 紊 为 2.5°x 
2.5°。 本 文 在 分 析 降 水 中 3D 和 8"”0 与 水 汽 来 源 的 关 
系 时 ,将 人 研究 期 划分 为 夏 半年 (5 一 10 月 ) 和 冬 半年 
(11 月 一 次 年 4 月 )”。 


ô, (7) 


稳定 同位 素 时 间 斥 度 效 应 及 水 汽 来 源 


2.3 研究 方法 
2.3.1 EK PSO 与 气温 和 降水 量 的 关系 ”气温 效 
应 指 降水 中 稳定 同位 素 组 成 与 气温 呈正 相关 关系 ; 
而 降水 量 效 应 指 随 着 降水 量 的 增 大 ,降水 中 稳定 同 
位 素 组 成 与 降水 量 呈 负 相 关 关 系 '“。 本 文 应 用 偏 
相关 系数 及 其 忆 值 ,研究 降水 中 8*0 与 气温 和 降水 
量 之 间 的 相互 关系 ”“”。 当 降水 中 8*0 的 气温 效应 和 
降水 量 效应 同时 存在 时 ,气温 效应 是 指 当 不 考虑 降 
水 量变 化 时 ,5"*0 随 气温 变化 的 一 阶 导数 ;而 降水 量 
效应 是 指 不 考虑 气温 变化 时 ,3*0 随 降水 量变 化 的 
一 阶 导数 ,并 符合 二 元 线性 模型 : 
8°O=a,t+a,7,+a,P,+é (8) 
SRB : a (fo °C!) All a (Do * mm |) Fy Be A EK 
量 效应 ;ao 为 常数 项 ;也 为 平均 气温 (% ) ;已 为 日 降水 
量 (mm); e 为 随机 误差 项 。 对 回归 方程 进行 方差 分 
析 和 检验 ,对 回归 系数 进行 :检验 ,判断 银川 平原 
是 否 有 显著 的 气温 或 降水 量 效 应 ,同时 应 用 回归 参 
数 w 和 的 偏 决定 系数 定量 气温 或 降水 量 效应 对 回 
归 方程 的 贡献 。 
2.3.2 降水 水 汽 来 源 后 向 轨迹 模拟 ”本文 利用 已 被 
广泛 应 用 于 水 汽 输送 等 方面 的 TrajStat 软件 ”" ,参考 
以 往 研 究 区 及 周边 地 区 水 汽 来 源 的 相关 人 研究， 
以 银川 站 为 轨迹 终点 ,以 6h 间 隔 ( 当 地 时 间 08:00、 
14:00、20:00 和 02:00) 模 拟 降 水 日 的 后 向 轨迹 ,轨迹 模 
拟 的 起 始 高 度 设置 为 1300 m、2500 m 和 5000 m3。 
由 于 平流 层 水 汽 平 均 停留 时 间 为 10 4 左右 2， 
此 将 回溯 时 间 设 定 为 240 h(10 d)。 因 本 文 只 关注 
水 汽 输 送 路 径 的 来 癌 ,采用 角度 距离 (Angle dis- 
tance) 算 法 聚 类 得 到 不 同时 段 下 主要 的 水 汽 输 送 
路 径 。 
2.3.3 潜在 源 贡 献 分 析 法 (PSCF) PSCF 是 一 种 基 
于 条 件 概率 函数 来 识别 潜在 污染 源 区 的 方法 。 
本 文 利用 TrajStat 软件 提供 的 PSCF ,将 水 汽 输 送 轨 
迹 经 过 的 区 域 划 分 为 0.5"x0.$" 的 水 平 网 格 (人 7) , 结 
合 水 汽 d-excess 日 数据 进行 分 析 ,并 将 夏 半年 水 汽 
d-excess 34) (4 0.8% Fil & 7B 4F 7k 4, d- excess 平均 
{H-0.1%0F} HNN BEL?” , PSCF 计算 公式 为 : 
PSCF, = (9) 


n.. 


式 中 :i 分别 为 经 度 和 纬度 ;nj 为 网 格 (i, 站 内 的 所 
有 轨迹 节点 数 ;mi 为 网 格 (i,)) 内 超过 平均 水 汽 d- 
excess PJE AJALE 17 4 BW; PSCF, APS (ij) Aa 
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平均 水 汽 d-excess 国 值 的 轨迹 出 现 的 概率 。 
PSCF 值 越 大 ,表明 网 格 内 水 汽 输 送 轨迹 对 应 的 
下 垫 面 蒸发 较 其 他 网 格 对 应 的 下 垫 面 蒸发 更 强 
烈 。 当 翅 较 小 时 ,PSCF 会 有 很 大 的 不 确定 性 ,因此 
以 某 一 网 格 改 小 于 研究 区 域内 各 网 格 平均 滞留 时 间 
的 3 倍 为 边界 ,引入 权重 机 数 Wi’: 
1.00, n;>80 
0.70, 20 <n; < 80 


i” 10.42, 10<n,< 20 wa 
0.05, n,< 10 
对 PSCF 进行 加 权 计算 : 
WPSCE, = W, x PSCF, (11) 


式 中 :WPSCF; 为 网 格 (i, 丫 的 WPSCF 值 ;WPSCF 值 
为 减少 不 确定 性 后 的 PSCF 值 。 

2.3.4 浓度 权重 轨迹 分 析 法 (CWT) PSCF 只 反映 
某 个 网 格 中 超过 平均 水 汽 d-excess 浆 值 的 轨迹 所 占 
的 比例 ,并 不 体现 其 对 目标 网 格 的 降水 贡献 。CWT 
则 可 以 确定 水 汽 的 潜在 蒸发 源 区 对 目标 网 格 的 相 
对 贡献 大 小 。 通 过 计算 每 个 网 格 的 水 汽 d-excess 
权重 值 定 量 给 出 经 过 该 网 格 的 水 汽 输送 轨迹 对 应 
的 平均 权重 水 汽 d-excess 值 ,以 反映 各 网 格 水 汽 d- 
excess 值 的 大 小 潜在 蒸发 源 区 和 局 地 再 循环 水 汽 
对 水 汽 输送 轨迹 中 d-excess 值 的 影响 状况 ,计算 公 
式 为 : 


N 
È Ci” Og 
Owr E — (12) 
Oik 
k=1 


式 中 :CWT; 为 网 格 (i,j) 上 的 平均 权重 水 汽 d-excess 
值 (%o"m?);k 为 气 团 轨 迹 数 目 ;N 为 轨迹 总 数 ;Ci 为 
轨迹 经 过 网 格 (i, 门 时 对 应 的 研究 站 点 水 汽 d- 
excess 值 的 大 小 (%o) ;a 为 轨迹 在 网 格 (i, 门 内 的 
停留 时 间 。PSCF 中 所 用 的 权重 函数 WW; 同样 适用 于 
CWT, 定 义 WCWT 如下: 

WCWT,=W,xCWT, (13) 
式 中 :WCWT; 为 网 格 (i, 站 的 权重 轨迹 水 汽 d-excess 
值 的 贡献 程度 ;WCWT 值 为 减少 不 确定 性 后 的 CWT 
1E (%o0+m™) o 


3 结果 与 分 析 


3.1 降水 中 稳定 同位 素 的 时 间 变 化 特征 
降水 中 5D . 80 值 的 变化 范围 介 于 -123.5%o~ 
22.2%o .-12.8%0~5.0%o ,降水 量 加 权 平 均值 分 别 为 


— 31.7%0 、-4.0%o , d- excess 值 的 变化 范围 介 于 
-37.6%o~17.0%o , 降水 量 加 权 平 均值 为 -0.1%o( 图 
2)。 降 水 中 稳定 同位 素 值 具有 明显 的 季节 变化 , 冬 
半年 降水 量 加 权 平 均 3*0 值 (-1.2%o) 偏 正 于 夏 半 年 
的 降水 量 加 权 平 均 3*0 值 (-4.1%o)。 与 降水 中 5"*0 
的 季节 变化 相同 , 冬 半 年 降水 量 加 权 平 均 d-excess 
值 (7.5%o) 明 显 偏 正 于 夏 半 年 的 降水 量 加 权 平 均 q- 
excess 值 (-0.5%o)。 在 夏 半 年 , 远 距 离 水 汽 输 送 使 
得 降水 中 稳定 同位 素 值 偏 负 ;而 在 冬 半 年 ,水 汽 
来 源 于 湿度 小 .蒸发 强 的 中 高 纬度 大 陆 性 气 团 , 沿 
途 降 水 较 少 ,降水 中 稳定 同位 素 值 偏 正 ”。 银 川 站 
夏 半 年 和 冬 半年 的 降水 量 加 权 平 均 3*0 值 (-3.0%o 
和 -0.8%o) 均 明显 偏 正 于 灵 武 站 (-5.1%o 和 -1.3%o ) 。 
与 降水 中 580 的 季节 变化 相反 , 灵 武 站 夏 半 年 和 冬 
半年 的 降水 量 加 权 平 均 d-excess 值 (4.5%o 和 10.9%o) 
均 显著 偏 正 于 银川 站 (-5.9%o 和 -1.0%o)。 
3.2 局 地 大 气 降 水 线 特征 

由 图 3 可 以 看 出 ,银川 平原 全 年 局 地 大 气 降水 
线 (LMWL) ,与 全 球 大 气 降 水 线 (GMWL) 相 比 , 其 和 斜 
率 和 截 忠 均 偏 小 ,原因 是 研究 区 地 处 西北 内 陆 干 旱 
区 ,雨滴 在 非 饱 和 大 气 降落 过 程 中 受到 比较 强烈 的 
二 次 蒸发 效应 的 影响 拓 ” ,降水 中 的 稳定 同位 素 发 
生 不 平衡 分 馏 , 导 致 重 同位 素 富 集 ,LMWL 的 斜率 偏 
小 。 冬 半年 LMWL 的 斜率 和 截 距 大 于 夏 半 年 ( 表 
2) ,表明 夏 半年 降水 受 云 下 二 次 莹 发 作用 的 影响 程 
度 更 大 。2 个 站 点 LMWL 的 斜率 和 截 距 差异 明显 ， 
与 GMWL 相 比 斜率 和 截 距 均 偏 小 ,可 能 是 云 下 二 次 
蒸发 效应 的 强 弱 程度 和 局 地 水 汽 再 循环 作用 不 
同 。 夏 半年 LMWL 的 斜率 和 截 距 较 为 接近 ,而 冬 半 
年 则 差异 较 大 ,主要 是 由 于 夏 半 年 银川 和 灵 武 站 点 
较为 干燥 的 大 气 环 境 和 较 强 的 蒸发 过 程 ; 冬 半年 ， 
水 汽 来 源 于 湿度 小 .蒸发 强 的 大 陆 性 气 团 ,沿途 降 
水 较 少 ,两 除 效应 较 弱 ,降水 中 氧 氧 稳定 同位 素 值 
偏 正 ,使 得 2 个 地 区 LMVWL 的 斜率 和 截 距 偏 大 。 
3.3 气温 效应 和 降水 量 效 应 

在 全 年 和 夏 半年 时 段 下 ,降水 中 8"0 与 平均 气 
温 均 呈 显著 的 正 相 关 关 系 ( 表 3)。 夏 半年 时 段 的 气 
温 效应 (0.258+0.037)%o .大 于 全 年 时 段 的 气温 
效应 (0.211+0.031)%o*%C1。 然 而 ,在 冬 半 年 时 段 , 降 
水 中 8*0 与 平均 气温 不 存在 显著 的 正 相 关 关 系 , 这 
可 能 与 研究 区 在 该 时 段 下 受 中 高 纬度 大 陆 性 气 团 
的 影响 有 关 。 与 以 往 对 西北 地 区 及 宁夏 降水 中 8*0 
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图 2 银川 平原 降水 中 8D ,350 All d-excess 的 逐日 变化 
Fig.2 Daily variations of 8D, ôO and d-excess in Yinchuan Plain 
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注 :LMWL 为 局 地 大 气 降水 线 ,GMWL 为 全 球 大 气 降 水 线 。 
图 3 银川 平原 局 地 大 气 降水 线 


Fig. 3 Local meteoric water line in Yinchuan Plain 


HPA mE TE AAFC 28 VO A WP”? ;在 全 年 和 夏 
半年 时 段 ,降水 中 8"0 与 降水 量 呈 正 相关 关系 ,降水 
量 效 应 不 显著 ,可 能 是 由 于 降水 中 5"0 与 平均 气温 
和 降水 量 之 间 相 互 影响 。 因 此 ,进一步 利用 偏 相关 
系数 ,人 研究 降水 中 8"0 与 平均 气温 和 降水 量 之 间 的 
相互 关系 ( 表 4)。 当 平均 气温 保持 不 变 时 ,3 个 时 段 


下 ,降水 中 50 与 降水 量 的 相关 系数 仍然 大 于 0; 而 
当 降 水 量 保持 不 变 时 , 除 冬 半年 以 外 ,全 年 和 夏 半 
年 时 段 降水 中 50 与 平均 气温 仍 旦 显著 的 正 相 关 关 
系 ,这 说 明 研 究 区 存在 显著 的 气温 效应 ,而 在 全 年 、 
夏 半 年 和 冬 半 年 时 段 下 ,降水 量 效 应 不 明显 。 这 可 
能 是 银川 平原 距 海 较 远 .大陆 性 气候 显著 , 受 西 风 
气流 影响 较 明显 ,而 来 自 西 太平 洋 及 印度 洋 的 暧 湿 
气流 受 高 山 阻挡 难以 到 达 , 季 风 对 研究 区 的 影响 较 
弱 , 使 得 该 区 域 的 降水 量 效 应 不 显著 %。 已 有 研究 
发 现 ,在 年 尺度 上 ,宁夏 及 银川 地 区 降水 量 效应 不 
MEOT ,但 次 降水 事件 尺度 上 部 分 区 域 夏 季 可 能 
存在 一 定 的 降水 量 效应 。 

3.4 水 汽 来 源 分 析 

3.4.1 后 向 轨迹 聚 类 分 析 在 夏 半年 (网 4a~c)， 
1500 mm 高 度 下 西方 路 径 中 距离 输送 的 水 汽 合计 占 
FE 79.60% , 短 距 离 路 径 输 送 的 局 地 蒸发 水 汽 占 比 
20.40% ,2500 mm 高度 下 的 水 汽 输送 路 径 与 1500 m 
度 相似 , 均 以 西方 路 径 占 主导 ,其 中 ,中 距离 水 汽 输 
送 路 径 占 比 85.33%, 较 1500 m 高 度 增加 7.20% , 短 
距离 路 径 输送 的 局 地 蒸发 水 汽 仅 占 14.67% , 而 
5000 m 高 度 下 均 为 西方 路 径 中 距离 输送 的 水 汽 ;在 


202201.00071v1 


chinaXiv 


FHEWY 


ChinaXiv 合 作 期 刊 


表 2 银川 平原 大 气 降水 线 斜 率 与 截 距 的 季节 变化 


Tab.2 Seasonal variations of slope and intercept of local meteoric water line in Yinchuan Plain 


= 银川 灵 武 银川 平原 
斜率 IRIE R 斜率 IRIE R 斜率 IE R 
夏 半年 5.38 一 8.45 0.87 5.52 -10.41 0.71 5.43 -9.71 0.85 
冬 半年 9.55 0.60 0.99 8.26 7.29 0.94 9.10 5.08 0.96 
全 年 6.72 -3.60 0.84 7.33 0.58 0.88 6.79 -2.79 0.84 


E:R REAR. FEl 


表 3 不 同时 段 银 川 平原 降水 中 6*0O 与 气温 \ 降 水 量 的 相关 关系 


Tab.3 Correlation of “O in precipitation with temperature and precipitation in Yinchuan Plain in different periods 


+ 


季节 T T(P) Py PT) 
夏 半年 0.310(0.045 ) 0.327(0.037) 0.158(0.318) 0.191(0.231) 
冬 半 年 0.848(0.069 ) 0.790(0.210) 0.588(0.297) 0.351(0.649) 
全 年 0.320(0.028) 0.319(0.031) 0.151(0.311) 0.149(0.323) 


注 :7(P,) 表 示 当 保持 降水 量 (P,) 不 变 时 ,降水 中 80 与 7 的 相关 系数 ;P.(7T) 表 示 当 温度 (7) 不 变 时 ,降水 中 5"*O 与 P, 的 相关 系数 。 对 于 7,P 值 是 


拒绝 假设 相关 系数 小 了 


F 零 的 概率 ;对 


F 书 ,, 已 值 是 指 拒绝 假设 相关 系数 大 了 


PAE ABER 


表 4 不 同时 段 银川 平原 降水 中 5 "0 与 气温 和 降水 量 的 二 元 线性 回归 


Tab. 4 Binary linear regression of 5"O in precipitation with temperature and precipitation in 


Yinchuan Plain in d 


ifferent periods 


季节 常数 项 气温 效应 气温 效应 偏 决定 系数 RMSE R F 
夏 半年 -11.402 0.258+0.037' 0.324 3.905 0.085 2.899 
冬 半年 -10.533+0.075 0.473+0.210 0.740 4.056 0.508 3.067 
全 年 -9.990 0.211+0.031° 0.315 4.052 0.082 3.061 


+ 


注 :RMSE 为 均 方 根 误差 ,为 方差 齐 性 检验 ;* 表 示 在 0.05 水 平 上 显著 相关 。 


冬 半年 (图 4d~f) ,1500 m 高 度 下 西方 路 径 中 距离 输 
送 的 水 汽 合计 占 比 92.73%, 而 北方 路 径 短 距离 输送 
AY) Jey HE KE TK IAM H 7.27% ,2500 m 高 度 下 西方 路 
径 长 距离 输送 的 水 汽 占 比 22.02% ,中 距离 输送 占 比 
73.13%, 相 比 1500 m 高 度 减少 21.14%, 短 距离 输送 
占 比 2.77% ,北方 路 径 占 2.08% ,5000 m 高 度 下 均 为 
西方 路 径 长 距离 输送 的 水 汽 ,水 汽 输 送 路 径 较 1500 
m 和 2500 m 更 为 集中 。 在 夏 半 年 和 冬 半 年 ,银川 站 
降水 的 水 汽 来 源 主 要 为 西风 水 汽 , 同 时 受 局 地 蒸发 
水 汽 的 影响 ,这 与 张 百 娟 等 ”的 研究 结果 一 致 ,其 
中 在 冬 半年 受 西风 水 汽 的 影响 更 显著 ,局 地 蒸发 水 
的 影响 相 比 夏 半年 有 所 减弱 。 

3.4.2 PSCF 分 析 在 夏 半 年 (图 5a~c),1500m 高 


究 区 西南 部 ,另外 ,东部 区 域 也 有 一 定 贡献 ;在 冬 半 
年 (图 5d~f) ,1500 m 高 度 下 大 于 0.9 的 高 值 区 域 范围 
较 大 ,延伸 路 径 较 远 ,主要 分 布 在 研究 区 周边 地 区 和 
西部 ,2500 m 高 度 下 WPSCF 大 于 0.9 的 高 值 区 域 面 
积 明显 缩小 , 且 向 南 迁 移 ,主要 分 布 在 研究 区 周围 及 
西部 区 域 ,5000 m 高 度 下 WPSCF 大 于 0.9 的 高 值 区 
域 面积 持续 缩小 旦 分 布 集中 ,主要 分 布 在 研究 区 周 
边 地 区 和 西部 , 相 比 1500 m 和 2500 m 高 度 ,西方 路 
径 区 域 水 汽 传输 对 研究 区 降水 的 影响 可 能 减弱 。 

3.4.3 CWT 分 析 ”在 夏 半年 (图 6a~c),1500 m RE 
下 WCWT 值 超过 6 o'm 的 区 域 主要 分 布 在 研究 区 
西部 ,研究 区 东部 也 有 一 定 贡献 ,2500m 高度 下 
WCWT 值 超过 6 %o*m ”的 区 域 主要 分 布 在 西南 部 ， 


度 下 WPSCF 大 于 0.8 的 高 值 区 域 主要 分 布 在 研究 区 
西部 ,2500 m 高 度 下 WPSCF 值 在 0.4~0.5 的 潜在 贡 
献 源 区 面积 明显 缩小 ,主要 分 布 在 研究 区 南部 ,此 
外 ,东南 部 潜在 源 区 的 贡献 率 也 较 高 ,而 5000 mm 高 
度 下 大 于 0.4 的 高 值 区 域 明显 增多 ,主要 分 布 在 人 研 


而 人 研究 区 西部 大 于 8 Ao’ m” 的 高 值 区 域 显 著 缩 小 ， 
5000 m 高 度 下 WCWT 值 大 于 3 %o*m 5 的 高 值 区 显 
著 增 大 且 分 布 较 集中 ,主要 分 布 在 研究 区 西部 ,其 
次 是 东部 ;在 冬 半 年 (图 6d~f) ,1500 m 高度 下 无 明 
显 高 值 区 域 ,2500 m 高 度 下 大 于 -3 hom? 的 区 域 显 
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半年 -1500 m 


稳定 同位 素 时 间 斥 度 效 应 及 水 汽 来 源 


(d) 冬 半年 -1500 m 


(f) 冬 半 年 -5000 m 
N 


人 


23.42% 


sa 7.13% 
~ 
21.11% 


9, = 
26.06% 19.51% 977% 


图 4 银川 平原 降水 水 汽 来 源 后 向 轨迹 聚 类 


Fig.4 Backward trajectories cluster of moisture Source in Yinchuan Plain 


著 增 大 ,主要 分 布 在 研究 区 周边 地 区 ,此 外 ,人 研究 区 
西南 部 也 有 一 定 贡献 ,而 5000m 高度 下 大 于 0.1 %o* 
nm? 区 域 显著 缩小 ,上 且 分 布 较 分 散 , 主要 分 布 在 研究 
区 西部 。 


4 结 论 
本 文通 过 研究 银川 平原 不 同时 段 下 降水 中 稳 


定 同 位 素 主 要 的 环境 效应 水汽 来 源 、 潜 在 蒸发 源 
区 及 其 贡献 程度 ,得 出 以 下 主要 结论 : 


(1) 银川 平原 降水 氧 氧 稳定 同位 素 组 成 存在 明 
显 的 季节 性 变化 , 冬 半年 降水 氧 氧 稳定 同位 素 组 成 
(-38.6%o+51.6%o 和 -4.5%o+5.2%o) 明 显 偏 正 于 夏 半 
年 (-40.9%o+17.7%o 和 -5.7%o+3.0%o ) ;不 同时 段 下 ， 
夏 半年 LMWL 的 斜率 、 截 距 (5.43、-9.71) 低 于 冬 半 
年 LMWL 的 斜率 、 截 距 (9.10、5.08) 和 银川 平原 全 年 
LMWL HIRERE AEE (6.79 .-2.79) 

(2) 在 全 年 、 夏 半年 和 冬 半 年 3 个 时 段 下 ,扣除 
降水 量 的 影响 ,降水 中 8*0 与 平均 气温 仍 呈 正 相关 
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注 :WPSCF 为 网 格 (i,)) 的 WPSCF 值 。 
图 5 银川 平原 夏 半 年 和 冬 半 年 WPSCF 分 布 


Fig.5 Distribution of WPSCF in Yinchuan Plain during the summer and winter half-year 


关系 ,气温 效应 分 别 为 (0.473+0.210)%o + Co! 
《0.258+0.037)%o.% 和 (0.211+0.031)%o ;扣除 
气温 的 影响 ,3 个 时 段 下 ,降水 中 30 与 降水 量 仍 呈 
正 相 关 关 系 ,说 明 降 水 量 对 降水 中 5*0 的 影响 不 明 
显 ,降水 量 效 应 不 显著 。 


(3) 夏 半 年 和 冬 半 年 ,人 研究 区 降水 受 西风 水 汽 
影响 显著 ,同时 受 局 地 蒸发 水 汽 的 影响 。 夏 半年 ， 
研究 区 周边 地 区 及 西部 西南 部 和 南部 是 主要 的 水 
汽 潜在 蒸发 源 区 ; 冬 半 年 ,研究 区 周边 地 区 和 西部 
是 主要 的 水 汽 潜在 蒸发 源 区 。 
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(a) 夏 半 年 -1500 
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(d) 冬 半年 -1500 m 
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轨迹 水 汽 贡 献 程度 。 


图 6 银川 平原 夏 半 年 和 冬 半年 WCWT 分 布 


Fig.6 Distribution of WCWT in Yinchuan Plain during the summer and winter half-year 
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Time scale effect of hydrogen and oxygen stable isotopes in precipitation 


and source of water vapor in Yinchuan Plain 


GAO Yang', HAN Lei™, HAN Yonggui’, HUANG Xiaoyu', PENG Ling*, LIU Lili’ 
(1. School of Geography and Planning, Ningxia University, Yinchuan 750021, Ningxia, China; 2. Joint International Research 
Laboratory of China-Arab Featured Resources and Environmental Governance in Arid Region, Yinchuan 750021, Ningxia, China; 
3. Ningxia Key Laboratory of Resource Assessment and Environment Regulation in Arid Region, Yinchuan 


750021, Ningxia, China; 4. School of Agriculture, Ningxia University, Yinchuan 750021, Ningxia, China) 


Abstract: The environmental effects of stable isotope and water vapor precipitation sources in the Yinchuan 
Plain, Ningxia Hui Autonomous Region, China, were studied using backward trajectory clustering analysis, the 
potential source contribution function, and the concentration- weighted trajectory analysis method to resolve vapor 
sources and potential evaporation source areas. Our results showed that (1) there were obvious seasonal variations 
in the hydrogen-oxygen stable isotope composition of precipitation in the Yinchuan Plain, and the hydrogen-oxygen 
stable isotope composition of precipitation in the winter half- year (—38.6%0+51.6%o0 and —4.5%0+5.2%0) was 
significantly more positive than that in the summer half-year (—40.9%0417.7%o0 and —5.7%0+3.0%o); The local 
meteoric water line was greatly improved. The slope and intercept of the LMWL in the summer half-year (5.43, 
—9.71) were lower than that in the winter half-year (9.10, 5.08), and that of the whole year (6.79, —2.79). (2) The 
temperature effect was significant at three time periods: (0.473+0.210) %o0-°C ', (0.258+0.037)%o-°C | and (0.211 
0.031)%o+C™' under the entire year, summer and winter half-year, respectively, the precipitation effect was not 
significant. (3) During the summer and winter half-years, the vapor sources of precipitation in the study area were 
primarily westerly water vapor, but were also influenced by local evaporative vapor. In the summer half-year, the 
evaporative vapor source areas were mainly distributed in the area surrounding the study area and the western, 
southwestern, and southern parts of the study area; During the winter half-year, they were mostly found in the 
area surrounding the study area and to the west. The study’ s findings may provide a theoretical foundation for 
selecting methods for identifying water vapor sources in arid regions using stable isotope techniques. 

Key words: atmospheric precipitation; stable isotope effect; water vapor source; potential evaporation source; 
different time scales; Yinchuan Plain 


